Magnetization plateaus in dimerized spin ladder arrays 
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We investigate tlie ground state magnetization plateaus appearing in spin-i two- leg ladders built 
up from dimerized antiferromagnetic Heisenberg chains and dimerized zig-zag interchain couplings. 
Using both Abelian bosonization and Lanczos methods we find that the system yields rather unusual 
plateaus and exhibits massive and massless phases for specific choices or "tuning" of exchange 
interactions. The relevance of this behavior in the study of NH4CUCI3 is discussed. 
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I. INTRODUCTION 

The study of low dimensional antiferromagnets is play- 
ing a major role and bringing new insights in our current 
understanding of collective spin fluctuations Part of 
the fascination of this field is that in the extreme quan- 
tum limit of 5* = 1/2 there is a rather complex cross-over 
between 1-d Bethe ground states and 2-d long range or- 
der. By now it is well established that, at least within 
the strong coupling regime, ladders with an even number 
of chains exhibit massive phases with purely short-range 
spin correlations. In contrast, odd-chain ladders display 
in general massless spin excitations resembling typical 
1-d ground states with power-law decaying correlations. 
A wealth of experimental investigations have confirmed 
these expectations for an increasing number of materials, 
such as families of low dimensional cuprates like Sr-Cu- 
O and La-Cu-0 compounds which can be well de- 
scribed in terms of spin-1/2 Heisenberg antiferromagnets 
confined to ladder geometries. 

Surprisingly, recent measurements of magnetization 
curves in NH4CUCI3 crystals at high magnetic fields ||] 
revealed rather unexpected features which contrast with 
the general behavior of spin excitations expected for 
these systems Specifically, this two- leg ladder 

S = 1/2 dimerized compound displays two magnetiza- 
tion plateaus at one-quarter and three-quarter of the sat- 
uration magnetization, irrespective of the external field 
direction. These results confront one of the central is- 
sues regarding the condition of fractional quantization for 
the appearance of massive spin excitations or magnetiza- 
tion plateaus under external varying fields. On general 
grounds this condition can be shown to be H|- [B| 



pN S{1~ {M)) (=, Z , 



(1) 



where p, N and S stand respectively for the periodic- 
ity of the ground state, the number of coupled chains 
and the total spin per site, whereas (M) denotes the 
system magnetization normalized to its saturation value. 
Thus, at sufficiently low temperatures two coupled spin- 
1/2 dimerized chains {p = N ^ 2) should exhibit, if 
any, two plateaus at (M) — and 1/2. However, ob- 
servations at temperatures down to 0.5°K and magnetic 
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fields up to 35 T carried out in |^ showed no evidence of 
these latter plateaus. The situation is rather intriguing 
since apart from Sd-effects in NH4CUCI3 as well as low- 
temperature structural transitions which can not be ruled 
out, the magnetization of a number of low dimensional 
halides and pnictides, e.g. KCuClg 1,10, CsaCuCU 
|]ll| , seems to be in fair agreement with Eq. (|l|). 

Motivated by these conflicting observations, recent 
theoretical studies [|l2|,^ pointed out that two-leg lad- 
der antiferromagnetic systems can exhibit vanishing spin 
gaps depending on the manner in which the array of cou- 
pling exchanges is realized. The key issue is that the 
interplay between massive generating mechanisms such 
as dimerization and interchain couplings eventually yield 
no cost in energy to create spin excitations under mag- 
netic fields. Therefore, the system switches from one 
ground state to another (short or long ranged), depend- 
ing in a highly non-trivial way on the particular choice 
or "tuning" of coupling exchanges. This is an ubiqui- 
tous sign of the importance of quantum fiuctuations of 
individual spins on the ground state. Depending on the 
exchange interactions, fluctuations can manifest them- 
selves collectively into many possible ground states, par- 
ticularly in lower dimensions where their effects are most 
pronounced. 

In this work we further analyze these findings by means 
of two independent and complementary techniques such 
as Abelian bosonization Q- and Lanczos diagonal- 
ization |jl^ of finite systems. We shall consider ladders 
of richer topologies including both frustrated interactions 
and dimerization along the interchain and intrachain cou- 
plings. Specifically, we focus attention on two dimerized 
spin-1/2 chains interacting through an isotropic Hamil- 
tonian of the form 
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where the Sn denote spin-1/2 operators associated to site 



The array of coupling exchanges are set as J, 
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J,j^j, and parameterized by Jn^ = Ji [1 +(— 
say for chain (1). To maintain purely antiferromagnetic 
interactions throughout the 2L spins of the ladder with 
periodic boundary conditions, the dimerization parame- 
ters are kept bounded by \5i\ < 1 ,i = 1,2. The corre- 
sponding zig-zag structure is schematized in Fig. 1. 

Besides its theoretical interest at room tempera- 
ture this triangular topology is actually realized both in 
NH4CUCI3 II and KCuClg |,§ crystals. As it was re- 
ferred to above, their magnetization curves exhibit quite 
different features thus, we are especially interested to in- 
clude an homogeneous field term of the type 
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so as to elucidate the combined effect of the above kinds 
of coupling arrays and applied magnetic fields h, say 
along the z direction. Also, it is worth pointing out that 
anisotropic exchanges can be included straightforwardly 
in our bosonization procedure as well as in the numerical 
analysis. 

In studying the ground state regimes of different ex- 
change parameter sets, two subcases of particular interest 
arise immediately |l^,|l|. Clearly, by setting S2 = —1 (1) 
we obtain a non frustrated dimerized system with stag- 
gered (plain) bond alternation as indicated in Fig. 2. The 
analyses given in Ref. [ p^|jT^ have shown that the mag- 
netization behavior of the former case resembles that of a 
single dimerized chain, as opposed to the non-staggered 
or plain situation {S2 — I ) , in which a net magnetization 
plateau shows up at (M) = 1/2 . It will turn out that our 
general zig-zag ladder interpolates continuously between 
the above scenarios and yields rather robust plateaus at 
(M) = and 1/2. However, the fine tuning of exchanges 
can suppress these massive regimes and furthermore, for 
particular subsets of the parameter space, in addition two 
magnetization plateaus at (M) =1/4 and 3/4 emerge si- 
multaneously. This latter issue actually could shed light 
on the measurements reported in Ref. 

The layout of the paper is organized as follows. In Sec- 
tion 2, we recast the low-energy excitations involved in 
each chain of Eq. (H) in terms of a conformal field theory 
of a free bosonic field compactified at a magnetization 
dependent radius. The dimerization along chains (i^i) 
and zig-zag interactions (^2) render the bosonization ap- 
proach particularly suitable to examine respectively weak 
{.J2/J1, < 1) and strong (J1/J2, \52\ < 1) coupling 
regimes. They are treated in turn in Section 2 A and 
2 B. Section 3 complements the magnetization behav- 
ior conjectured by the analytic approach in a variety of 
non-perturbative scenarios. An exact numerical treat- 
ment of magnetization contours for finite systems up to 
2L = 24 spins is given using a recursion-type Lanczos 
algorithm applied on each ma gne tization subspace. 
Standard extrapolation procedures |15| to the thermody- 
namic limit then enable an independent test of the results 
obtained via bosonization techniques. We end the paper 



with Section 4 which contains our conclusions, along with 
some remarks on pros and limitations of the present work. 



II. ABELIAN BOSONIZATION 

Following a recent analysis discussed as in Refs. 0- [D, 
we will apply the by now standard method of Abelian 
bosonization at Hamiltonian (j^). In this formalism an 
antiferromagnetic homogeneous chain is described by a 
compactified free bosonic field (Z)*^"-* whose dynamics is 
governed by 

Hi-) dx {yK{dJi-)f + ^{d.,4-)f) (4) 

The dual field is defined as usual 11 = 8x4'^°') ■ v is 
the Fermi velocity and the Luttinger constant K, which is 
a function of the magnetization (M^"') and an eventual 
XXZ anisotropy A, governs the conformal dimensions 
of the bosonic vertex operators and can be obtained ex- 
actly from the Bethe ansatz solution of the XXZ chain 
(see e.g. |^ for a detailed summary). It is related to the 
compactification radius i? of |^ by K^^ = 2nR^. 
In terms of these fields, the spin operators read 
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a : cos(2fc>x -I- 72^0^"^) : + 

g±,ia) ^ (—1)^ ■ g±»V2ir0<"' 

(bcos{2k''p^x + V2^(t>'^'') 



(5) 
(6) 



where the colons denote normal ordering with respect to 
the ground state with magnetization (Af*^"^). The Fermi 

momentum fc^'' is related to the magnetization of the 

a*'' chain as k^p'^ = (1 - (M('')))7r/2. The effect of an 
XXZ anisotropy A and/or the external magnetic field 
is then to modify the scaling dimensions of the phys- 
ical fields through K = K{(Mi'''>) , A). The magnetiza- 
tion also modifies the commensurability properties of the 
spin operators through fc^?, as can be seen from (H), (|). 
The non-universal constants a, b and c can be in general 
computed numerically (see e.g. |l^, for the case of zero 
magnetic field) and in particular the constant c has been 
obtained exactly in 

Notice that the inclusion of an XXZ anisotropy in our 
study is not only motivated to generalize the analysis but 
primarily by the fact that, for non-zero magnetization, 
the SU{2) symmetry is broken from the beginning. As we 
shall see in Section 2 A, the particular SU{2) symmetric 
case, (A = 1, (A/) = 0), has to be analyzed differently 
since our analysis breaks explicitly this symmetry. We 
address the reader to Ref. [Q where the full symmetric 
was analyzed using a formulation in terms of Majorana 
fermions. 
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A. Weak interchain coupling regime 

Here we take a = J2/J1 <C 1 and (5i ^ 1. In this 
regime |(52| <C 1 corresponds to a weakly coupled two-leg 
zig-zag ladder made up of dimerized chains, whereas for 
1 the system approaches a staggered (52 ~ —1) or 
plain (^2 = 1) ladder dimerized as in Fig. 2(a) and Fig. 
2 (b) respectively. 

Using Eqs. (||, ||), the low energy Hamiltonian can be 
written in this regime as 



X 

A2 ^ cos \^kFX -f \/27r(0(^) 4 

X 

+ A3^cos[\/2^(,/)(i) 

X 

+ A4^cos [\/2^(0(i) - 

X 

+ As ^(-l)'^+"cos [2/cF(a; + l/2) 

a;, (a) 

V2^0('')(a;)j 



+ Ae X! X!'^'^^ 2fcFa: -t- \/2^0^''^(x) 

(a) = l a: 



(7) 



where 



A2 , j = 1 , 

A2 [(1 - <52) + (1 + ^2) cos(2fc;^)] , J = 2, 3 , 
Aj/a oc <( (52, j = 4, (8) 
^1 , j = 5 , 
I (M),J = 6. 

In the last two expressions the proportionality factors are 
non-universal functions of the XX Z anisotropy A2, the 
magnetization (A/) and bi. 

In the above equation we have suppressed marginal 
terms, due to the presence of more relevant interactions. 
Also in the case of non-dimerized zig-zag interactions, 
there are parity breaking terms, discovered in |l8| , which 
should be analyzed differently since they have non-zero 
conformal spin. However, by explicitly including these 
latter terms in the RG computations, it can be shown 
that they do not change the conclusions presented here. 
The S't/ (2) symmetric case (A = 1, (M) = 0) has been 
studied in for the case of zero zig-zag dimerization. 

Let us analyze the structure of the magnetization 
curves predicted by this effective Hamiltonian. The rel- 
ative field 0_ = 01 — 02 is always massive due to the 
fact that the relevant perturbation terms A3^4 are al- 
ways commensurate. The analysis for the diagonal field 
0+ = 01+02 is more subtle so we will consider separately 
different values for the magnetization. The possible val- 
ues of the magnetization at which plateaus can appear 
are given by the general expression (|]) jJJ^. However, 



due to the presence of frustration in certain region of 
the parameter space, possibly Eq. (|l|) should include 
an extra factor 2 to account for a possible enhancement 
in the periodicity of the ground state. This restricts 
the set of possible magnetization plateaus to the values 
(A/) = 0, 1/4, 1/2,3/4 (apart from saturation). 

For (M) = 0, (i.e. fc^ = 7r/2), the perturbation A2 
becomes commensurate and hence opens a gap for the 
diagonal field 0+. One could however argue that taking 
into account other (radiative) corrections (coming mainly 
from A5 and Ag) one could close the 0+ gap, simply by 
making the amplitude of this perturbation term to van- 
ish. This naive analysis has been confirmed by numeri- 
cal computations jl^ , showing that there exists a whole 
curve in the parameter space where this indeed happens. 
In particular, for ^2 = — 1 the critical line turns out to 
be 

J2 



oc 



'1 • 



(9) 



This phenomenon was originally suggested in |]l2| using 
non-linear sigma model techniques and has been the also 
studied numerically in . 

For (M) = 1/2, (ie. kp ^ 7r/4), the plateau can 
open due to the radiatively generated (relevant) terms, 
(coming from A5 and Ag), which are of the form 

o?h{M)j{h,8^, A2)(-l)" cos(4A:j.a; + V2^0+), (10) 

where /((5i, (^2, A2) vanishes for the case of (^2 — — 1 , 
as already pointed out in ||l3|. As it will be shown in 
Section 3, possibly there is another point at which the 
(Af) = 1/2 plateau closes. The existence of such effects 
can be predicted using the bosonization formalism, but 
the precise location of this point, however, cannot be 
obtained due to the presence of non-universal constants 
in Eqs. d), d). 

The case of (Af) = 1/4,3/4, (i.e. kp = 37r/8,7r/8), is 
less clear since the commensurate operators that can be 
generated to open this plateau are irrelevant. One possi- 
ble candidate is the operator (—1)"^ cosi^kpx -I- 2\/27r0+) 
which is irrelevant. One can speculate whether this oper- 
ator can become relevant in some region of the parameter 
space by taking into account the Ai perturbation which 
changes the compactification radius of the 0+ field, and 
hence the dimensions of the perturbation terms. How- 
ever, in first approximation the above mentioned opera- 
tor would become relevant for values of the couplings far 
from the region where our approach can be considered 
valid. 



B. Strong interchain coupling regime 

Let us consider now the opposite regime, i.e. a 3> 1 . 
In this regime it is convenient to rewrite the two-leg zig- 
zag ladder Hamiltonian as that of a single chain with al- 
ternating intrachain coupling J2(l±(52) and next-nearest- 
neighbor (NNN) interactions Ji(l ± 8\) but alternating 
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every two sites. Namely, NNN spins at (4n, 4n + 2) and 
(4n + 1, 4n + 3) pair locations are coupled by Ji (1 — Ji) , 
whereas spins at {An — 2,4n) and (4n — 1 , 4n + 1 ) interact 
through Ji(l + 6i) . This is illustrated in Fig. 3. 

Using the same approach as described above, in the 
present case we get an effective Hamiltonian for a single 
bosonic field perturbed by the following terms 

X 

H- 72 (^) COS |^2fe_FX + \/27r(/)(x) 

X 

X 

+ ^74(x)cos Akpx + 2V2Tr (l){x) . (11) 

X 

Here 71,2,4(2:) are proportional to Si and 73 to 52- The 
terms proportional to Si have an extra modulation due to 
the alternation of the NNN couplings referred to above, 
which hereafter we call "two-by-two" modulation. 

For 61 = and (M) — only the last term is com- 
mensurate and it has dimension 2 in the isotropic case 
(A = 1). This term is the one responsible for the opening 
of the gap at zero magnetization in the NNN antiferro- 
magnetic Heisenberg chain, where, as it is well known, 
the gap opens at a critical value of the NNN coupling 
through a Kosterlitz Thouless (KT) transition. Its effect 
would have been the same in the present case, but due 
to the presence of the dimerization Si along the chain, a 
gap will always be present. This will be corroborated in 
Section 3. 

For (A/) — 1/2, the perturbation 72 survives the con- 
tinuum limit due to the extra "two-by-two" alternating 
factor, and will hence be responsible for the plateau at 
this value of the magnetization. 

The case of (M) = 1/4, 3/4 is again more subtle in this 
limit, and it can be seen that the operator that could 
be responsible for the appearance of these plateaus is 
generated from a combined effect of the chain dimeriza- 
tion ((—1)^(52) and the two-by-two alternating pari of the 
NNN exchange (±5i Ji). The operator generated through 
this mechanism is proportional to cos(2\/27r 0(x)), which 
is irrelevant, and could become relevant at certain critical 
value of the coupling (cx JiSiS2)- Again, this KT point 
is not reachable within our perturbative approach. 



III. NUMERICAL ANALYSIS OF FINITE 
SYSTEMS 

To enable an independent check of the magnetization 
scenario obtained within the bosonization approach, we 
now turn to a numerical finite-size analysis of the original 
ladder Hamiltonian (0). A number of numerical studies 
of triangular ladders (or equivalently, of Heisenberg chain 



with NN exchanges), have been reported already ( [p[|22|] 
and references therein). However, the effect of dimer- 
ization along both interchain and intrachain couplings, 
which is crucial for the appearance of non-trivial magne- 
tization plateaus, yet requires further numerical efforts. 

We focus attention on the ground state energy ob- 
tained from an exact diagonalization of finite systems via 
a recursion type Lanczos algorithm [Q applied on each 
magnetization subspace with — |0, 1, ...,L}. Since 
the magnetic field considered in Eq. (^) is coupled to the 
conservation of S^/L = (M), we can readily relate the 
energy per spin Ch at finite fields to those at /i = just by 
taking eh = eg — h (M). Thus, all results addressed be- 
low were obtained from computations with h = . Also, 
to avoid unwanted effects introduced by both the ladder 
topology and periodic boundary conditions, even multi- 
ples of ladder lengths up to i — 12 spins were taken 
throughout. 

The huge dimensionality of the spaces involved, grow- 
ing as {^^), constrained us to use in the heaviest situa- 
tions, i.e. 5'^ = 0, 1 , at most seven Lanczos vectors per 
tridiagonalization cycle ||l^. Nevertheless, the numeri- 
cal accuracy was kept bounded by 10~''/i/Ji employing 
typically up to forty cycles of recursion. On the other 
hand, the rather small number of Lanczos vectors used 
in the computations allowed for an efficient management 
of a complete, vector by vector reorthogonalization. As 
is known jlj], this latter procedure becomes crucial to 
avoid the emergence of spurious eigenvalues caused by 
machine rounding errors which tend to build up expo- 
nentially with the number of iterations, no matter what 
precision is used. In what follows we limit our analysis 
to isotropic coupling exchanges, though preliminary cal- 
culations including anisotropy in the field direction yield 
qualitatively similar results. 

We begin by examining the validity of Eq. (|l|) and test 
the "fine tuning" effects conjectured within the bosoniza- 
tion analysis. Upon setting Si = 0.7, S2 = 0.475 and 
J2/J1 — 1, the application of the Lanczos procedure to 
the triangular ladder Hamiltonian (|^) yields a quite un- 
usual behavior. Specifically, in addition to the typical 
(Af) = and 1/2 plateaus, two massive phases appear 
simultaneously at (M) = 1/4 and 3/4 as is shown in 
Fig. 4. Rather unexpectedly, at this point of the pa- 
rameter space the spin ladder seems to decouple into a 
quasi-four level system (see the region near the vertical 
lines of Fig. 5). In principle this renders size effects 
less pronounced, possibly due to some enhanced or hid- 
den symmetry whose origin is yet difficult to elucidate. 
Nevertheless, our results provide a strong numerical ev- 
idence indicating that zig-zag structures can yield non 
trivial magnetization plateaus such as those observed in 
NH4CUCI3 crystals §. 

Interesting as it is, we are actually motivated to obtain 
the whole magnetization curve of this latter compound. 
Therefore, we turn to the scanning of the exchange pa- 
rameter space for different regimes. Figs. 5(a), 5(b) and 
5(c) display respectively typical magnetization contours 
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or "phase diagrams" for weak (J2/J1 = 0.5, 5i = 0.9 ), in- 
termediate (J2/J1 = l,6i = 0.7), and strong (J2/J1 — 5, 
di — 0.5), coupling regimes. This is a compact form 
of representing conventional magnetization curves for a 
wide range of interchain couplings. Here each line is asso- 
ciated to successive values of (M) which increase mono- 
tonically with the applied field h. For example, the mag- 
netization plateaus of Fig. 4 are contained completely 
within the vertical line of Fig. 5(b). In general, we found 
that the (M) = and 1/2 plateaus remain robust in a 
variety of scenarios, though their widths can be eventu- 
ally "fine tuned" to yield massless gaps. 

In studying the mass gap extrapolation towards their 
thermodynamic limits (i.e. the energy gap to create an 
excitation of total spin 5*=! asL^cx)),we fitted the 
whole set of finite-size results for 4 < i < 12 {L even), 
using a variety of standard procedures. These range from 
linear to logarithmic and van den Broeck-Schwartz type 
methodologies of convergence |l^ , which basically yield 
analogous results with at least two significant digits. We 
draw the reader's attention to panels 6(a)-(b), 6(c)-(d), 
and 6(e)-(f) in which we display respectively gap extrap- 
olations around (M) = 1/2 and corresponding to the 
ground state regimes exhibited in Figs. 5(a)-(c). Similar 
gap extrapolations for (Af) = 1/4 and 3/4 beyond the 
"symmetry" or finite size collapse region denoted by the 
vertical lines of Figs. 5(a) and 5(b), would be unreliable 
given the scarcity of available data. In our case, this is 
translated in the availability of matching sizes, namely 
L = 4, 8, 12, already constrained by the studied values 
of (M) . Thus, it remains unclear whether or not empty 
wide "bands" or plateau regions for (M) = 1/4 and 3/4 
could be actually present in Figs. 5(a)-(c). 

To complement the analysis of finite (vanishing) (M) = 
1/2 gaps for ^2 = 1 (-1) given in Section 2 A, we see 
that their widths remain stable upon setting 62 < ^ , 
{62 > — 1 )• Moreover, the wide minima of Figs. 6 (a), 
(c) , and (e) suggest an infinitely continuous ground state 
transition at 62 = —1, which is in line with the KT sin- 
gularity conjectured in Section 2 B. In fact, our data 
strongly support this picture for a variety of coupling 
regimes, at least within the region 62 G (—0.9 , —0.25) 
where finite size effects are less pronounced. This can 
be observed from the semi-log representation of the data 
displayed in Fig. 7. Thus finally, the triangular ladder 
turns out to interpolate smoothly between the staggered 
{S2 = —1) and plain (62 — 1) dimerization arrays of the 
square ladders studied in ||l3| . 

IV. CONCLUSIONS 

In this work we have studied the magnetization phase 
diagram of a two leg zig-zag ladder with dimerization 
both along the legs and the zig-zag coupling, by means 
of Abelian bosonization methods complemented by Lanc- 
zos diagonalization of finite clusters up to 24 spins. From 



the bosonization analysis we conclude that the (M) = 
plateau is robust and it is present in the full parameter 
space, except on a certain zero-measure set. We have 
confirmed numerically these expectations and obtained 
the above massless excitations as displayed in the lower- 
most panels 6(d) and (f). 

Spin excitations around (A/) = 1/2, also turn out to 
be massive and robust as they show up in the whole pa- 
rameter space, except in the limit (52 — — 1 and for a 
certain curve J2/J1 = f{Si,S2)- This latter feature is 
observed in Figs. 6(a), (c) for weak and intermediate cou- 
pling regimes near ^2 = and around S2 ^ 0.5 for strong 
coupling regions as shown in Fig. 6(e). 

Regarding the issue of (M) = 1/4,3/4 plateaus, they 
are observed only within a fine-tuned region, possibly 
bearing an enhanced symmetry of the Hamiltonian (see 
e.g. Figs. 4 and 5). Their appearance is hard to predict 
using Abelian bosonization techniques. 

Finally, though there arc intermediate values of S2 ca- 
pable of closing gaps around both (M) = and 1/2, after 
scanning a representative set of the parameter space, we 
could not find any evidence of common closing points so 
as to explain the suppression of these plateaus at a time 
in NII4CUCI3 crystals. In principle, this would rule out 
the ladder Hamiltonian (^) as a suitable model to account 
completely for the experiments reported in [|| . Nonethe- 
less, we trust that its success to describe the simultaneous 
emergence of rather unusual plateaus at (M) = 1/4 and 
3/4 (at least near the data collapse region of Fig. 5), will 
make our system worth to consider as an antecedent for 
future studies in that direction. 
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Ji(l + 6i) 




FIG. 1. Schematic view of doubly alternating zig-zag lad- 
ders showing both interchain (Ji), and intrachain (J2) ex- 
change couplings along with their respective dimerization pa- 
rameters (5i and 82 ■ 
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FIG. 2. Representation of non-frustrated ladders with 
dimerized chain bonds Ji(l ± 5i) = J(l ± 5) and interchain 
coupling J' = 2J2, obtained by setting (a) ^2 = — 1 (stag- 
gered dimerization) and, (b) 62 = ! (plain dimerization). 
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FIG. 3. Chain representation of dimerized zig-zag ladder 
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FIG. 5. Critical fields of dimerized zig-zag ladders for (a) 
J2/J1 = 0.5 with 5i = 0.9, (b) J2/J1 = 1 with Si = 0.7, 
and, (c) J2/J1 = 5 with (5i = 0.5. Solid, dashed and short 
dashed lines denote respectively the results for L = 12,8 and 
4, whereas J — max {Ji, J2} . Vertical linos in (a) and (b) 
indicate regions where 1/4 and 3/4 plateaus emerge at a time. 
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FIG. 6. Width of magnetization plateaus for (M) = 1/2 

in (a), (c), (e) and, (M) = in (b), (d), (f). Data of first, 
second and tiiird panel refer respectively to J2/J1 = 0.5 with 
3i = 0.9, J2/J1 = 1 with Si = 0.7 and, J2/J1 = 5 with 
Si = 0.5 . Solid lines in descending order indicate respectively 
results for L = 6, 8, 10 and 12, whereas lowermost dashed 
curves [slightly visible in (b) and (f ) ] denote gap extrapola- 
tion to the thermodynamic limit. 
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FIG. 7. Gap width extrapolations for (A/) ~ 1/2 suggest- 
ing a KT singularity oc e~^^^^~^^^ on approaching S2 = — 1. 
Triangles, squares and circles denote respectively weak, in- 
termediate and strong coupling regimes displayed in turn in 
Figs. 6 (a), 6(c) and 6(e). Straight lines are guides to the eye. 
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